Direct photons at forward rapidities in high-energy pp collisions 
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We investigate direct photon production in pp collisions at the energies of RHIC, CDF and LHC, at different 
rapidities employing various color-dipole models. The cross section peaks at forward rapidities due to the 
abelian dynamics of photon radiation. This opens new opportunities for measurement of direct photons at 
forward rapidities, where the background from radiative hadronic decays is strongly suppressed. Our model 
calculations show that photon production is sensitive to the gluon saturation effects, and strongly depends on 
the value of the anomalous dimension. 

PACS numbers: 13.85.QK,13.60.Hb,13.85.Lg 



I. INTRODUCTION 

Photons radiated in hadronic collisions not via hadronic de- 
cays are usually called direct. They carry important informa- 
tion about the collision dynamics, not disturbed by final state 
interactions. In particular, the hadronization stage is absent, 
so the theoretical interpretation is simpler that in the case of 
hadron production. A unified description for radiation of vir- 



tual (Drell-Yan reaction) and real photons within the color 
dipole approach was proposed in (ll, 0]. This description 
does not need to be corrected either for higher order effects 
(K-factor, large primordial transverse momentum), or for the 
main higher twist terms'. The corresponding phenomenol- 
ogy is based on the universal dipole cross section fitted to DIS 
data and provides a rather good description of data, both the 
absolute normalization and the transverse momentum depen- 
dence |4]. Predictions of the inclusive direct photon spectra 
for the LHC at midrapidity and the azimuthal asymmetry of 
producedprompt photons in the same framework are given 
in Refs. IS la]. Comparison with the predictions of other ap- 
proaches at the LHC can be found in Refs. 

Intensive study of the dynamics of hadronic interactions 
at high energies and search for signatures of nonlinear QCD 
effects, like saturation |9], or color glass condensate (CGC) 
flOil , have led to considerable experimental progress towards 
reaching smallest Bjorken x. The typical experimental set up 
at modern colliders allows to detect particles produced in hard 
reactions in the central rapidity region, while the most ener- 
getic ones produced at forward (backward) rapidities escape 
detection. The first dedicated measurements of hadron pro- 
duction at forward rapidities, by the BRAHMS experiment in 
deuteron-gold collisions at RHIC disclosed an interest- 
ing effect of nuclear suppression, which can be interpreted 



' Some higher twist corrections, which are specific for forward rapidities, 
still have to be added j^. 



as a breakdown of QCD factorization II12I1 . The zero degree 
calorimeters detecting neutral particles, neutrons and photons, 
at maximal rapidities are being employed at RHIC and are 
planned to be installed at LHC. One could increase the trans- 
verse momentum coverage of these detectors by moving them 
away from the beam axis. Through this note we would like to 
encourage the experimentalists to look up at these opportuni- 
ties. 



The important advantage of measurements of direct pho- 
tons at forward rapidities is a significant enhancement of the 
signal-to-background ratio. Indeed, the photons radiated by 
the electric current of the projectile quarks, which stay in the 
fragmentation region of the beam, form a bump at forward 
rapidities (see Fig. 5). At the same time, gluons are radi- 
ated via nonabelian mechanisms by the color current which 
flows across the whole rapidity interval. Therefore gluons 
are radiated mainly in the central rapidity region ifisll . and are 
strongly suppressed in the beam fragmentation region. Such a 
suppression is even more pronounced for hadrons from gluon 
fragmentation, and for photons from radiative decays of those 
hadrons. Another source of background photons, hadrons pro- 
duced via fragmentation of the valence quarks is also strongly 
suppressed due to the shift towards small fractional momenta 
related to the quark fragmentation function and the kinematics 
of radiative decay. Thus, direct photon production at forward 
rapidities should be substantially cleared up. 



Here we perform calculations for direct photon production 
at large pj and various rapidities in proton-proton coflisions 
at the energies of RHIC, Tevatron and LHC. We employ the 
color dipole approach and compare the predictions of several 
contemporary models for the dipole cross section, based on 
the idea of gluon saturation. 
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II. PHOTON RADIATION IN THE COLOR DIPOLE 
FORMALISM 

Production of direct photons in the target rest frame should 
be treated as electromagnetic bremsstrahlung by a quark in- 
teracting with the target, as is illustrated in Fig.[T]in the single 
gluon approximation, which should be accurate at large trans- 
verse momenta of the photon. Since the quark trajectories 
before and after photon radiation have different impact pa- 
rameters, and the corresponding terms in the bremsstrahlung 
amplitude have different signs, one arrives at an expression, 
which is formally identical to the amplitude of an inelastic 
dipole-target interaction 1 1]. This is only a formal procedure 
of calculation, while no real dipole is involved in the process 
of radiation. 

Calculation of the transverse momentum distribution is 
more involved since the direct and complex conjugated 
amplitudes correspond to incoming quarks with different im- 
pact parameters. It substantially simplifies after integration 
over transverse momentum of the recoil quark, so one is left 
with two dipoles of different sizes ?i and r^, and the cross 
section gets the factorized form Hi, 



d{lna) d^pT 



X e 



iPT-iri-h) 



(|)* (a, r 1 )(|)y^ (a, F2 )Ey(x2 , n , r2 , a) , ( 1 ) 



where n and ri are the quark-photon transverse separations 
in the direct and complex conjugated amplitudes respectively; 
a = Py / denotes the fractional light-cone (LC) momen- 
tum of the radiated photon. Correspondingly, the transverse 
displacements of the recoil quarks in the two amplitudes are 
ari and ar2 respectively. The LC distribution amplitude for 
the qj Fock component with transverse separation r has the 
form. 



2k 



(2) 



Here X;,/ are the spinors of the initial and final quarks and 
Ko{x) is the modified Bessel function. The operator O has the 
form, 

d^imfa?e*-{nxa) + ae* • (a x V) ~ /(2 - a)e* • V , (3) 

where e is the polarization vector of the photon; n is a unit 
vector along the projectile momentum; and V acts on rj. The 
effective quark mass serves as infra-red cutoff parameter, 
which we fix at niq « 0.14GeV, since all dipole parametriza- 
tions considered in this paper also assume the light quark mass 
equal to 0.14 GeV. 

In equation ([T]) the effective dipole cross-section Ly is a lin- 
ear combination of the qq dipole-proton cross sections, 

1 



Ey(x2,n,r2,a) = -{Ogg{x2,ari)+Ogg{x2,ar2)} 



Ogg{x2,a{ri-r2)). 



(4) 
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FIG. 1 : Direct photon radiation by a quark interacting with a gluon 
in the target. 



Here and in what follows xi 2 are the Bjorken variable of the 
beam and target partons. 



(5) 



and r| is the photon rapidity in the cm. of pp collision. 

Since only quarks and antiquarks can radiate photons, the 
hadronic cross section is given by the convolution of the par- 
tonic cross se ction Eq. ^ with the proton structure function 



^2 



dG{pp yX) 
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daFl' 



da{qp yX) 
dad^PT ' 



dxfd^ Pt x\ +X2 

(6) 

where Xf — xi — X2 is the Feynman variable. This relation 
needs commenting. The transverse momentum distribution 
of quark bremsstrahlung should be convoluted with the pri- 
mordial transverse motion of the projectile quark. Differently 
from the parton model, in the dipole approach one should rely 
on the quark distribution function taken at a soft scale. Evolu- 
tion to the hard scale is performed via gluon radiation, which 
is encoded into the phenomenological dipole cross section fit- 
ted to DIS data for the proton structure function. Since the 
quark primordial motion with a small (soft) mean transverse 
momentum does not affect the photons radiated with large pr 
||4|], we neglect the transverse momentum convolution and use 
the integrated quark distribution. 

However, a word of caution is in order The dipole cross 
section Ogg{x2,r) includes gluon radiation which performs the 
evolution and leads to an increase of the transverse mo- 
mentum of the projectile quark. However, it misses the 
evolution of the xi -distribution, which i s esp ecially important 
at forward rapidities (see discussion in llql ). since the quark 
distribution falls off at xi — ^ 1 much steeper at high Q^. In 
order to account for this effect and provide the correct xi- 
distribution, we take the integrated quark distribution in ^ 
at the hard scale. 

We use the parametrization for the proton structure function 
given in Ref. iITtIi . Following Ref. |4] , for the hard scale 
entering the proton structure function in Eq. we choose 

Q = PT- 
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FIG. 2: The total dipole-proton cross section aqq(x,r) at a fixed 
X = 10^^ in various dipole models introduced in Sec. III. 



III. MODELS FOR THE DIPOLE CROSS SECTION 

The dipole cross-section is theoretically unknown and 
should be fitted to data. Several parametrizations proposed 
in the literature are employed here to investigate the uncer- 
tainties and differences among various models. 

A popular parametrization proposed by Golec-Biernat and 
Wusthoff (GBW) model fl?] is based on the idea of gluon 
saturation. This model is able to describe DIS data with the 
dipole cross-section parametrized as, 



oo 1 



(7) 



where jc-dependence of the saturation scale is given by 
= {xq/x)^ GeV^. The parameters Oo = 20.1 mb, xq — 
5.16 X 10^"*, and A, = 0.289 were determined from a fit to F2 
without charm quarks. A salient feature of the model is that 
for decreasing x, the dipole cross section saturates at smaller 
dipole sizes. The saturation scale in the GBW reduces with 
the inclusion of the charm quark ifTgl] . After inclusion of the 
charm quark with mass = 1.5 GeV, the parameters of the 
GBW model changed to Oo = 23.9 mb, jcq = 1 ■ 1 1 x 10"'* , and 
A, = 0.287. Both parametrization sets give a good description 
of DIS data at x< 0.01 and [0.25,45] H. 

One of the obvious shortcoming of the GBW model is that it 
does not match the QCD evolution (DGLAP) at large values of 
Q^. This failure can be clearly seen in the energy dependence 
of aJof for > 20 GeV^, where the the model predictions 
are below the data lfT8ll20|] . 

A modification of the GWB dipole parametrization model, 
Eq. (|7]i, was proposed in Ref. |20] by Bartels, Golec-Biernat 



and Kowalski (GBW-DGLAP) 



GBW-DGLAP 



{x,r) —<3q[ 1 — exp 



% r a.s {/J )xg{x,ij ) 
3oo 



where the scale /j^ is related to the dipole size by 



2 c 



(8) 



(9) 



Here the gluon density g{x,/j^) is evolved to the scale /j^ 
with the leading order (LO) DGLAP equation. Moreover, 
the quark contribution to the gluon density is neglected in the 
small X limit. The initial gluon density is taken at the scale 
Ql = IGeV^ in the form 



xg{x,i? 



— 



(1 



\5.6 



(10) 



where the parameters C ~ 0.26, <3q = 23 mb,/jQ = 0.52 GeV^, 
Ag = 1.20 and = 0.28 are fixed from a fit to the DIS data 
forx<0.01 and in a range of ^ [0.1,500] GeV^ 120]. The 
dipole size determines the evolution scale /j^ through Eq. 
The evolution of the gluon density is performed numerically 
for every dipole size r during the integration of Eq. ([T]). There- 
fore, the DGLAP equation is now coupled to our master equa- 
tion ([Hi. It is important to stress that the GBW-DGLAP model 
preserves the successes of the GBW model at low and its 
saturation property for large dipole sizes, while incorporating 
evolution of the gluon density by modifying the small-r be- 
haviour of the dipole size. 

Since the linear DGLAP evolution may not be appropriate 
for the saturation regime, lancu, Itakura and Munier proposed 
an alternative color glass condensate (CGC) model [21,1 , based 
on the Balitsky-Kovchegov (BK) equation ll22ll . The dipole 
cross section is parametrized as. 



rQs < 2 



: rQ, > 2 

(11) 

where Q, = Qs{x) = {xo/xf^l^ GeV, Y = ln(l/x), and K = 
X"{'Ys)/x'{'Ys) where % is the LO BFKL characteristic func- 
tion. The coefficients A and B in the second line of ( fTTl i are 
determined uniquely from the condition that a(x,r), and its 
derivative with respect to rQs, are continuous at rQs = 2: 



<7? 



B 



(1-%) 



(12) 

The parameters Ys = 0.63 and K = 9.9 are fixed at the LO 
BFKL values. The others parameters iVo ~ 0.7, Oo = 35.7 mb, 
xo = 2.7 X 10"^ and A. = 0. 177 were fitted to F2 for x < 0.01 
and Q^ < 45 GeV^ and including a charm quark with = 1 .4 
GeV. Notice that for small rQs < 2, the effective anomalous 
dimension 1 — js in the exponent in the upper line of Eq. ( fTTI) 
rises from the LO BFKL value towards the DGLAP value. 

It should be stressed that this CGC model is built based 
on the solution of Ref. ll23ll for r Qs > 2 and a form of the 
solution for rOs < 1, but in the vicinity of r 0= 1 /g.s it is given 
in Refs. ll24ll25ll. 
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FIG. 3: Inclusive direct photon spectra calculated with various color- 
dipole models at mid-rapidity T| = at the energies of RHIC, ^/s = 
200 GeV (upper panel), and Tevatron, -Js =1.8 TeV (lower panel). 
The NLO QCD curve is from the authors of reference (30ll (given 
in table 3 of Ref. fsill ). Experimental data are from the PHENIX 
experiment fs?] at r] = 0, and from the CDF experiment (sTl [33ll 
at |t|| < 0.9. The error bars are the linear sum of the statistical and 
systematic uncertainties. 



Notice that calculation of the /^j- -distribution given by 
Eq. ([TJ needs only knowledge of the total dipole cross sec- 
tion and is independent of the impact parameter dependence 
of the partial elastic dipole-proton amplitude. Nevertheless, 
we consider also the model proposed by Watt and Kowalski 
ll26ll . Although the main focus of this model is the impact 
parameter dependence (b-CGC), which is irrelevant for our 
calculations, the integrated cross section is different from the 



,fb-CGC 



{x,r)=2 



(13) 



where a^^'^{x,r,b) is given by Eq. (fTTl i with the saturation 
scale Qs which now depends on impact parameter. 



^ r 
1 



exp 



2BcGc/ 



1 

2W 



(14) 



The parameter Bqg = 7.5GeV^ is fitted to the f-dependence 
of exclusive J/^ photoproduction. It has been shown that 
if one allows the parameter js to vary along side with other 
parameters (in contrast with CGC fitting procedure where js is 
fixed with LO BFKL value), it results in a significantly better 
description of data for F2 with the value of js — 0.46, which is 
remarkably close to the value of js ~ 0.44 recently obtained 
from the BK equation IztIi . Other parameters obtained from 
the fit are: 5Vd = 0.558, xq = 1.84x 10"*' and A. = 0.119. 

In order to demonstrate the importance of saturation, we 
will also use a non-saturated model (No Sat) fitted to F2 with 
X < 0.01 and G [0.25,45] GeV^: 



(15) 



where Qs is defined in Eq. (fT4l i. The parameter Jeff is de- 
fined for rgs < 2 as Jeff = + 1^ 7§7' ™d for rQ, > 2 
as Jeff = Js- The other parameters are given by js = 0.43, 
5Vd = 0.568, xo = 1.34 X lO"*^ and X = 0.109 |26]. Surpris- 
ingly, the fit obtained with such an oversimplified model is as 
good as the other models with x^/d.o.f. = 0.92 (although it 
should certainly fail to explain data on diffractive DIS Il28fl . 
which are sensitive mainly to large size dipoles). Notice that 
we use here the No Sat model for a qualitative argumentation. 

At first sight this result could be used as an argument that 
the data is not sensitive to the saturation effect. However, the 
actual meaning of this exercise is quite opposite. It is well 
known that the saturation effects start being essential when the 
anomalous dimension reaches the value jcr = 1 ^ Jeff ~ 0.37 
(see Refs. H M, Si)- We will show that at very forward 
rapidities at LHC, the diffusion term in the anomalous dimen- 
sion is less important. Therefore, what we actually demon- 
strate is that the value of the anomalous dimension should be 
larger than jcr at very forward rapidities for LHC, and because 
of this the saturation effects have to be taken into account. 

The second comment on this model (see Ref. |26]) is that 
it is actually a model which contains saturation, and the dif- 
ference with the CGC model (see Eq. ( fTTT )) is only one: this 
model is written for dipoles with sizes close to 1 /Q^. In- 
deed, comparing Eq. (fTTT i and Eq. ( fTSI ) one can see that they 
treat differently the region rQ^ > 1. The CGC model de- 
scribes this region as solution to the BK equation deeply in 
the saturation region |23], with a phenomenological matching 
at rQs = 2, while this model uses the solution to the BK equa- 
tion ill Hi, m for rQs > 2 but r close to l/gj. Therefore, 
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FIG. 4: Direct photon spectra obtained from various dipole model at 
midrapidity (upper panel) and forward rapidities (lower panel) at the 
LHC energies for pp collisions. 



it is not appropriate to call this model "no saturation model", 
nevertheless, we use this name as a terminology. 



Summarizing, we can claim that direct photon production is 
sensitive to saturation effects. In conclusion, the success of the 
so-called 'no saturation model' can be interpreted such that at 
the LHC we will be still sensitive to the kinematic region close 
to the saturation scale. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In Fig. [3] we compare predictions of various dipole models 
with data for inclusive prompt-photon production from RHIC 
at = 200 GeV and from the Tevatron at = 1 .8 TeV. A 
word of caution is in order here. All the above parametriza- 
tions for the dipole cross-section have been fitted to DIS data 
at X < 0.01. This coiTesponds to pr < 2 GeV at the RHIC 
energy , so the PHENIX data plotted in the upper panel of 
Fig.[3]are not suited for a model test. Notice that the CDF data 
plotted in the bottom panel of Fig. [3] were obtained with a so 
called isolation cut, which is aimed at suppression of the over- 
whelming background of secondary photons originated from 
radiative hadron decays. This might change_the cross-section 
within 10 - 20% for the CDF kinematics (Q. 

One can see from Fig. [3] that various dipole models pre- 
sented in the previous section with explicit saturation give 
rather similar results at small pr- At high pr, CGC, b-CGC 
and GBW-DGLAP models which incorporate QCD evolution 
provides a better description of data compared to the GBW 
model. The No Sat model with the diffusion term, defined in 
Eq. ( fTSl ), gives similar results as the b-CGC model for both 
RHIC and CDF energies (not shown in the plot). In order to 
understand the role of the diffusion term in the the anomalous 
dimension, we show in Fig. [3] the results with two extreme 
limits Ye// = 0.43, 1. It is clear that the dipole model without 
explicit saturation as given by Eq. ( fTSI ) with y^-// = 0.43, does 
not describe the data either from PHENIX, or from CDF (we 
do not show in Fig. [3] No Sat. with Jeff — 0.43 curve for CDF, 
since it is about two orders of magnitude above the other mod- 
els and data). However, changing the anomalous dimension to 
the DGLAP value with Jeff = 1, dramatically changes the re- 
sults and brings the curves (No Sat model) at both energies 
of RHIC and Tevatron (at r| = 0) within the ranges of other 
dipole models with saturation and the No Sat model in the 
presence of the diffusion term. Therefore, the diffusion term 
in the anomalous dimension is very important at both RHIC 
and CDF energies. 

We have recently shown that the color dipole formulation 
coupled to the DGLAP evolution provides a better description 
of data at large transverse momentum compared to the GBW 
dipole model |4]. In Fig.lH upper panel, we show the predic- 
tions of the GBW (with charm quark) and the GBW-DGLAP 
models for LHC energies ^/s — 5.5,14 GeV at midrapidity 
for the transverse momentum nptopr— 200 GeV. In Fig. H] 
lower panel, we show the predictions of various color-dipole 
models for ^/s —14 GeV at different rapidities. Generally, the 
discrepancy among predictions of various models at moderate 
Pt is not very large. This can be also seen from Fig. |5] where 
we compare, as an example, the GBW and b-CGC models, 
which expose very different structures (see Eqs. ( I7I13I ) and 
Fig. (|2|). In the inserted plot in Fig. |4] we show the effect 
of unitarization within the GBW model, namely using the ex- 
ponent in Eq. (|7| as a dipole cross-section (r2 model). One 
can also see that the discrepancy between the GBW and the 
r2 model increases at forward rapidities, though it is still not 
appreciable. 

In Fig. m we also show the results for the model without 
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FIG. 5: The difference between the GBW and the b-CGC models 
for the LHC energy -/s = 14 TeV at different rapidities. Inserted 
plot: the discrepancy between the GBW and the r2 model (the GBW 
model without explicit saturation) at various rapidities for =14 
TeV in pp collisions. 



saturation. At the midrapidity again the results of the No Sat 
model with the DGLAP anomalous dimension 1 — y^ff = 
is close to other saturated models. At the same time, at very 
forward rapidities, the anomalous dimension 1 — Jeff = 0.57 
which is close to the value predicted from the BK equation 
iIztIi . will be in favour of other models. At forward rapidi- 
ties, the diffusion term in the anomalous dimension is not 
important more, since it gives similar result as with a fixed 
Jeff = 0.43. This indicates that direct photons production at 
different rapidities at LHC is rather sensitive to the saturation. 
Again, since the values of anomalous dimension turn out to be 
larger than jcr — 0.37, such a description of the experimental 
data indicates at a large saturation effect. 

In Fig.|6] the differential cross-section of photon radiation 
at the energy of LHC is plotted versus rapidity at fixed trans- 
verse momenta pr — 1,3,5 GeV. Calculations were performed 
with several models for the dipole cross section. All of them 
lead to a substantial enhancement of the photon production 
rate at forward rapidities. One can see that the larger the sat- 
uration scale is, the stronger is the peak. From Fig. |6] it is 
again obvious that No Sat model at very forward rapidity is 
within the dipole model's predictions with an explicit satu- 
ration. However, the peak disappeared since at about mid- 
rapidity Jeff = 0.43 is too small. In principle, the peak can be 
also reproduced in the No Sat model if one allows an anoma- 
lous dimension running with energy and transverse momen- 
tum pt- The appearance of the peak at forward rapidity is a 
direct consequence of the abelian nature of the electromag- 
netic interaction. In the case of gluon radiation, the peaks 
at forward-backward rapidities is be replaced with a kind of 
plateau at central rapidities, which is indeed observed in data 
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FIG. 6: Invariant cross-section for direct photon production in pp 
collisions as a function of rapidity calculated with various color 
dipole models r| for various fixed pj shown in the plots. 



for hadron production. 

To conclude, by this letter we encourage measurements of 
direct photons at forward rapidities in pp collisions at mod- 
ern colliders . These experiments will be a sensitive tool for 
search for saturation effects, since they will allow to access 
the smallest possible values of Bjorken x in the target. Be- 
sides, the background of photons from radiative hadronic de- 
cays should be significantly suppressed. As we demonstrate 
in Fig. |6] direct photons are enhanced, even form a bump, at 
forward rapidities. At the same time, gluon nonabelian ra- 
diation is known to be strongly suppressed in this region, so 
hadron and decay photons are also suppressed. We provided 
predictions for the cross section of direct photon production 
at various rapidities for pp collisions at LHC employing dif- 
ferent models for the dipole-proton total cross section. 
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